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Abstract 

By using first-principles calculations, we propose that ZrSiO can be looked as a three-dimensional 
(3D) oxide weak topological insulator (TI) and its single layer is a long-sought-after 2D oxide TI 
with a band gap up to 30 meV. Calculated phonon spectrum of the single layer ZrSiO indicates it 
is dynamically stable and the experimental achievements in growing oxides with atomic precision 
ensure that it can be readily synthesized. This will lead to novel devices based on TIs, the so called 
“topotronic” devices, operating under room-temperature and stable when exposed in the air. Thus, 
a new field of “topotronics” will arise. Another intriguing thing is this oxide 2D TI has the similar 
crystal structure as the well-known iron-pnictide superconductor LiFeAs. This brings great promise 
in realizing the combination of superconductor and TI, paving the way to various extraordinary 
quantum phenomena, such as topological superconductor and Majorana modes. We further find 
that there are many other isostructural compounds hosting the similar electronic structure and 
forming a WHM-iaiaSiy with W being Zr, Hf or La, H being group IV or group V element, and 
M being group VI one. 
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I. INTRODUCTION 


Topological insulator (TI) has nontrivial band topology characterized by Z 2 topologi¬ 
cal invariant under time-reversal symmetry (TRS).^®^ Similar to the systems with integer 
quantum Hall effect, the bulk band structure of TIs is insulating with unavoidable edge or 
surface states on the boundary connecting the valence and conduction bands. Especially for 
the two-dimensional (2D) TI, it can host quantum spin Hall effect (QSHE) with ID helical 
edge states, where the states moving in opposite direction have opposite spin. Therefore, 
the backscattering is prohibited as long as the scattering potential does not break TRS and 
such helical edge states provide a new mechanism to realize non-dissipative electronic trans¬ 
portation, which promises potential application in low-power and multi-functional device 
based on TI. From this point of view, a 2D TI with large band gap, chemically stable under 
ambient condition exposed in the air, easy to prepare and consisting of cheap and nontoxic 
elements is highly desired.®^ The study of TI based quantum devices would foster a field 
of “topotronics”. Such 2D TI is also important and crucial to the realization of topological 
superconductivity and Majorana modes through proximity effect.^ 

However, up to now there are only a few systems, which have been proved to be 2D 
TiEEIsuch as the quantum well of HgTe/CdT^^ and InAs/GaSb.^ Both of above systems 
require very sophisticated and expensive MBE growth and only show QSHE at ultra-low 
temperature. These extreme requirements have obstructed further studies and possible 
applications of 2D TIs. The theoretical efforts in Ending and designing 2D TI candidates 
keep going on and some of them are listed here: silicend^, chemically decorated single layer 
honeycomb lattice of Sn,^G^and Bi or Sb,^i3HIEl^ buckled Bil^, single layer ZrTes, HfTe^ 
and Bi 4 Br 4 ,^ transition-metal dichalcogenide (TMD) in IT® and haeckelitd^D® structure. 
But none of them has been confirmed experimentally. Recently, Weng et al proposed a large 
band gap 2D TI in oxygen functionalized MXene,® which brings the hope in realizing TI in 
oxygen contained compounds, which are naturally antioxidant and stable upon exposed in 
the air. This also stimulated the effort in this paper in searching for new oxide TIs. 

In this paper, by using first-principles calculations we demonstrate that the oxide material 
ZrSiO is a weak TI and a single layer of it is a 2D TI with achievable band gap up to 30 
j 44 eY[ 23 EZl have calculated its edge states and analyzed the physical mechanism for band 
inversion. We have found many other isostructural compounds with formula of WHM 
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(W=Zr, Hf or La, H=Si, Ge, Sn or Sb and M=0, S, Se and Te) possessing the similar 
electronic strnctnre. 


II. COMPUTATIONAL DETAILS 

We have employed the Vienna ab initio simnlation package (VASPp^*^ for most of the 
density fnnctional theory (DFT) based first-principles calcnlations. Exchange-correlation 
potential is treated within the generalized gradient approximation (GGA) of Perdew-Bnrke- 
Ernzerhof type® Spin-orbit conpling (SOG) is taken into acconnt self-consistently. The 
cut-off energy for plane wave expansion is 500 eV and the k-point sampling grid in the self- 
consistent process was 12x12x3 for 3D case and 12x12x1 for the single layered structure. 
The crystal structures have been fully relaxed until the residual forces on each atom is 
less than 0.001 eV/A. A vacuum around 16 A between layers was used to minimize the 
interactions between the layer and its periodic images. The possible underestimation of 
band gap within GGA is checked by non-local Heyd-Scuseria-Ernzerhof (HSE06) hybrid 
functionapi® calculation. To explore the edge states, we build a tight-binding model for a 
slab around 120 unit-cell thick with edges along a or 6 lattice.® The maximally localized 
Wannier functions (MLWF)®® of d orbitals of W and p orbitals of H are generated as the 
basis set by using the software package OpenMX.®® 


III. RESULTS AND DISCUSSION 

Crystal Structure. The bulk ZrSiO belongs to a big family of ZrSiS- (or PbFGl-) type 
compoundP^ with chemical formula of WHM (hF=transition metal or rare earth element; 
H, M=main group elements). There are over 200 members in this family,® and the crys¬ 
tal structure in Fig. [T] can be looked as stacking of hve square nets of hF, H and M in 
sequence of [... M-W-H-W-M ...] along c axis. The relatively weak M-W bonding con¬ 
necting two neighboring slabs show in Fig. leads to quasi layered structure of them. In 
fact, it was thought that the cross over from 3D to 2D in ZrSiM (M=0, S, Se and Te) 
happens between M=Se and Te.® Nevertheless, this gives out a natural cleavage plane 
for getting 2D layer of these compounds and we will show the band dispersion along c is 
quite narrow than that of in-plane for several members in this family. It is noticed that 
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LiFeAs ZrSiO 


FIG. 1: (Color online) Crystal structure of (a) LiFeAs (b) ZrSiO. (c) Top view of single 
layer ZrSiO and 3D Brillouin zone with high symmetrical crystal momenta indicated. The 
weak bonds between Si-0 have been indicated by dashed lines, which form a natural 
cleavage plane to get 2D single layer ZrSiO. 


the isostructural pnictide NbSiAs has been found to be superconductor.®^ The intensively 
studied iron-pnictide superconductor LiFeAs also has the iso-structure but with the cation 
and anion sites being switched to form the so called anti-PbFCl (or anti-ZrSiS) type struc- 
ture.®^The space group of ZrSiO is P4/nmm (No. 129). The theoretically optimized lattice 
constants are a=6=3.316 A and c=7.354 A, being underestimated by about 6-7% compared 
with experimental values a=6=3.52 A and c=7.93 A.®' The is within the common error 
in lattice constant estimation from DFT calculation within GCA. The optimized Wyckoff 
site of atoms are Zr (0.25, 0.25, 0.308), Si (0.75, 0.25, 0) and O (0.25, 0.25, 0.6112), also 
consistent with the experimental values. The general features of the band structure from 
the optimized crystal structure and the experimental one are quite similar. To be consistent 
in the description of the main text, the theoretically optimized crystal structure is used in 
the following discussion. 

Band structure of 3D ZrSiO. The band structure of 3D ZrSiO is shown in Fig. [^a). 
There are two distinct features. One is the apparent band crossing points around the Fermi 
level indicating the existence of band inversion as found in many other materiald®^ with 
nontrivial topological band structures.® The second feature is its band dispersion in 
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plane (F-X-M-F) resembles that in kz='K plane (Z-R-A-Z) and there is no band inversion 
happen along F-Z direction. This means that its electronic structnre is essentially more 2D 
like than 3D. The further group theory analysis indicates that the bands crossing in both 
kz='i^ and 0 planes are protected by glide mirror perpendicular to the z-axis because the 
two bands forming the crossing points belong to different eigenvalues of that glide mirror. 
As discussed in Ref. H21 - H51 these band crossing points form a closed node-line circle in 
the /c^=0 and A;^=7r plane, respectively, when SOC is not considered. Including SOC into 
calculation opens a gap along the node-lines turning 3D ZrSiO into a system with hnite 
gap at all k points but compensated electron and hole pockets at the Fermi level. As we 
have discussed in Ref. 031 and 011 the crossing points on the node-line are not exactly on the 
same energy level due to the particle-hole asymmetry. There will be a hnite global band gap 
instead of compensated electron-hole pockets. Since the material has inversion center, we 
can determine its Z 2 invariances through the production of the parities at the time reversal 
invariant momenta (TRIM), which gives (0; 001) indicating that it is a weak TI if the band 
gap opening on the node-lines due to SOC is larger enough to open the indirect gap. ESM] 
The above results also indicate that 3D ZrSiO can be viewed as the stacking of 2D TI along 

c-axis,^2SI25] ZrTe^ and Bi4Br4.®' 

2D TI of Single layer ZrSiO. As discussed above, we take a single layer of ZrSiO by 
breaking the weak Si-0 bonds. In order to get the binding energy between layers of ZrSiO, we 
calculate the total energies of 3D and single layer ZrSiO with the same setting of parameters 
in calculations. The binding energy is about 1.166 or 1.503 eV as shown in Table |T] and 


III calculated within GGA and LDA, respectively. It is much bigger than the interlayer 


binding energy of graphite and ZrTe^ and Bi 4 Br 4 .^. In this point of view, the interlayer 
interaction in ZrSiO is weak covalent bonding rather than van der Waals interaction. 

On the other hand, the band structure of single layer ZrSiO essentially has the same band 
topology as that in any plane perpendicular to kz axis in 3D ZrSiO. The only difference 
is in the strength of particle-hole asymmetry. On considering these and the experimental 
technique of well controlled layer-by-layer growth of oxided^, we think the single layer ZrSiO 
is readily synthesized. We perform full geometrical optimization of the free standing single 
layer ZrSiO. The cell parameters a and b have been relaxed to a=6=3.225 A. In the direction 
c, the atoms Zr move close to Si layer by about 0.0251 A, and atoms O move close to Si 
layer by about 0.004 A. Such small relaxation can also be ascribed to its quasi-2D feature. 
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FIG. 2: (Color online) Band structure for (a) bulk ZrSiO within GGA and single layer 
ZrSiO within (b) GGA and (d) GGA+SOC. The inset in (b) indicates the Node-Line due 
to band inversion around Fermi level protected by in-plane mirror symmetry. The inset in 
(c) indicates the gap-opening in Node-Line when SOC is further included. 



FIG. 3: (Color online) The phonon spectrum of optimized single layer ZrSiO. 


The phonon spectrum of the fully relaxed structure is shown in Fig. It has no imaginary 
frequency and indicates the dynamical stability of single layer ZrSiO. 

The band structures shown in Fig. (b)-(d) are for fully optimized 2D ZrSiO. As shown 
in Fig. 1^ (b), the two Dirac points on the paths F-X and F-M are nearly at the same energy 
level close to Fermi energy. When SOC is further considered, a global band gap around 0.01 
eV is obtained as indicated by the total density of states (DOS) in Fig. (c) and the band 
structure in Fig. (d). The insets in Fig. (b) and (d) clear show the gap opening by SOC 
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FIG. 4: (Color online) The edge state of single layer ZrSiO along a or 6 lattice vector, (a) 
The Dirac cone like edge state connects the bulk valence and conduction bands, (b) The 
wave function at the sampled k in (a) on the edge state shows its decay into bulk as 

function of number of unit cell. 


along the node-line circle. We hnd that by slightly compressing the thickness of the slab 
or equivalently, expanding the in-plane lattice constants would decrease the particle-hole 
asymmetry and change the size of band gap. For example, taking the experimental in-plane 
lattice constant and compressing the ZrSiO layer by 4% would lead to a global band gap 
about 30 meV. 

Since it is central symmetric, the Z 2 number 1 is determined by counting the parity of all 
occupied states at four time reversal invariant momenta.l^^® This conhrms that the single 
slab of ZrSiO is a 2D TI. Its edge state along a or 6 lattice vector is shown in Fig. which 
is calculated using the Green’s function methocP^ based on the tight-binding hamiltonian 
in basis of Wannier functions of Zr d and Si p orbitals.^^siMl 'jpg edge states are characteristic 
Dirac cone like. We have also studied the decaying of edge state into the interior of 2D 
ZrSiO. As shown in Fig. (b), it decays quite slow compared with the surface states of 3D 
TIs of Bi 2 Se 3 or Bi 2 Te 3 . Therefore, a quite thick slab is needed, more than 80-unit-cell thick 
along a or 6 lattice, to get gapless Dirac cone along its edge. 

Band Inversion Mechanism. The crystal field splits five Zr d orbitals into d^'^-y'i, dxy, dz 2 , 
and dyz -l- d^z manifolds and three Si p orbitals into pz and Px + Py The DOS plot and the 
fatted bands in Fig. [^a) indicate that the bands involved into band inversion around Fermi 
level are mainly composed of Zr dx^-y-^ and Si Px + Py orbitals. In each unit cell, there are 
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FIG. 5: (Color online) The band inversion mechanism in single layer ZrSiO. (a) The fatted 
band structure for Zr dx 2 _y 2 and Si + Py orbitals, (b) The fatted band structure for La 
5dx2_y2 and Sb + Py orbitals, (c) and (d) are the band inversion mechanism with and 
without hybridization between dx' 2 -y 2 and px + Py orbitals, corresponding to ZrSiO and 

LaSbTe, respectively. 


two nonequivalent Zr and Si atoms. The band inversion mechanism can be understood as 
following. Firstly, the degenerate atomic orbitals form the bonding and anti-bonding states, 
with the splitting determined by the effective hopping parameters of td and tp. The bonding 
and anti-bonding states from dx 2 -y 2 orbitals constitute the 62 ^ (odd parity) and Big (even 
parity) representation of the little group at F. Those from degenerate Px + Py orbitals 
construct the 2D representation and Eg of It is noted that the dxz + dyz of Zr 

atoms form the same 2D representation as px -l- Py orbitals and they are heavily hybridized. 
Secondly, the banding effect along F-X leading to band broadening and band inversion. As 
we note that both the glide mirror {M^|(I, |, 0 )} (with ab plane as mirror plane) and the 
inversion symmetry P anti-commute with the screw rotation {( 72 x 1 ( 1 , 0 , 0 )} at three TRIM 
points at the zone boundary X, Y and M. As a consequence, similar to the situation in 
ZrTes we can prove that all the states at those TRIM points at the zone boundary are 
doubly degenerate with opposite glide mirror eigenvalue and parity, which leads to a unique 
property that the Z 2 invariance here is determined by the order of states at F point point, 
because all the contributions from the parity of the states at the zone boundary TRIM points 
always cancels each other. Since the states at X must be paired with glide mirror eigenstates 







































with opposite eigenvalues as shown in Fig. (c) and (d), when the bands evolves 

from r to X, there are unavoidable band crossings for two typical order of states at F from 
two different such materials (ZrSiO and LaSbTe). As discussed in the previous section, these 
band crossings are protected by the glide mirror symmetry and form the node-line inside 
of the mirror plane as shown in Fig. (b). With SOC included, the band gap will open in 
ZrSiO and its Z 2 invariance can be determined by the states at F point to be odd. The htted 
band structure from the tight-binding model with basis of and px -l- Py orbitals can 

well reproduce the one from first-principles calculation as shown in Fig. (c). The relative 
energy level of dx^-y2, Px + Py and Fermi level Ej determines the band topology. With the 
same number of electrons as ZrSiO, it is trivial insulator if dx 2 -y 2 bands are fully occupied as 
shown in Fig. (c). However, if dx 2 -y 2 are higher than px +Py orbitals, Fermi level will pass 
through the band crossing point and it becomes a 2D TI with SOC included. We find that 
the isostructural LaSbT^^ with the same number of electrons exactly follows this picture 
and its band structure is shown in Fig. (b). The La 5d orbitals are full empty and have 
one less valence electron than Zr or Hf, while the Sb has one more valence electron than 
Si and other group IV element. The single layer LaSbTe is a 2D TI with global band gap 
around 0.011 eV with SOC included. The band inversion and the band gap of ZrSiO and 
LaSbTe are further checked by HSE06 calculations in Fig. We hnd that in both cases the 
band inversion is well kept. The particle-hole asymmetry in ZrSiO is modihed a little and 
the global band gap is missing. For LaSbTe, the band gap is slightly increased to 0.015 eV. 

Single layer ZrSiO on substrate SrTiO^. In fact, the single layer material should be 
supported by substrate in experiment. To check the influence of substrate on the electronic 
structure of ZrSiO, we put a single layer of ZrSiO on commonly used substrate SrTiOa. 
There are two types of top surface for SrTiOs, namely, the SrO terminated surface and 
the Ti022 terminated one, as shown in Fig. (a) and (c), respectively. Three unit cells 
of SrTiOs are taken as substrate. When ZrTiO is put on each surface, the cations are put 
on top of anions to minimize the total energy and the a, b lattice constants are expanded 
to fit the SrTiOs substrate. The initial distance of them is about half of cell parameter of 
SrTiOs. There is a vacuum layer about 13 A is used. The three-layer atoms of SrTiOs at the 
interface and all atoms of single layer ZrSiO are fully relaxed. The geometrical structure has 
nearly no change. The fatted bands with weight on single layer of ZrSiO are shown Fig. 
(b) and (d). We hnd that the most features of free standing ZrSiO is well kept, especially 
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FIG. 6: (Color online) The band structures of (a) ZrSiO without SOC, (b) ZrSiO with 
SOC, (c) LaSbTe without SOC and (d) LaSbTe with SOC, calculated within GCA (black 

lines) and HSE06 (red lines). 


the band crossings around Fermi level are well reproduced inside of the band gap of SrTiOs 
substrate. On considering the underestimation of SrTiOs band gap by about 2.0 eV within 
GCA, we believe that SrTiOa is quite suitable for supporting ZrSiO. 


The WHM family. We have further calculated the other isostructural compounds WHM 
with IF=Zr, Hf, if=Si, Ge, Sn and M=0, S, Se and Te. The band structures for both the 
3D and single layer 2D WHM have been shown in Figs. i i ID inj in and All 
of these bands show the similar band structure around the Fermi level. The difference 
is in the strength of particle-hole asymmetry, which brings most of them are semimetals 
with overlapped electron-hole pockets. As we have discussed above, the expansion in a, 
b lattice constant and compression of the slab thickness can be helpful in decreasing the 
asymmetry and thus leads to the global band gap. From these hgures, we hnd the particle- 
hole asymmetry increases as M changes from O to Te. At the same time, the binding energy 
between layers with GGA in table [I] and [IT] decreases when M changes from O to Te. It 


is also validated by the example of ZrSiM with LDA in Table, Therefore, the further 
detailed works on the WHM family are necessary and left for successive paper. 
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FIG. 7: (Color online) The geometrical structure and band structure of single layer ZrSiO 
on (a)-(b) SrO terminated surface and (c)-(d) Ti02 terminated surface of SrTiOs. The red 
circled lines are the fatted bands with the circle size proportional to the weight on single 

layer of ZrSiO. 


IV. CONCLUSIONS 

Based on the hrst-principles calculations, we have predicted that 3D ZrSiO is a weak 
TI and its single layer ZrSiO is an oxide 2D TI with band gap around 30 meV. The other 
members in this family noted as WHM with hF=Zr, Hf, i7=Si, Ge, Sn and M=0, S, Se 
and Te have the similar band structure and are most probably 2D TIs. The isoelectronic 
compound LaSbTe has been shown share the similar band structure and share the same 
underlying physics in band inversion. 
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FIG. 8: (Color online) The band structure for (left) 3D and (right) 2D ZrSiM with M=0, 
S, Se and Te. The experimental crystal structure of them are used. 
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FIG. 9: (Color online) The band structure for (left) 3D and (right) 2D ZrGeM with 
M=0, S, Se and Te. Except M=0 case, the others are calculated with experimental 
structure. The experimental crystal structure of ZrSiO is used for ZrGeO. 
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FIG. 10: (Color online) The band structure for (left) 3D and (right) 2D ZrSnM with 
M=0, S, Se and Te. The experimental crystal structure of ZrSiO, ZrGeS and ZrGeSe are 
used from ZrSnO, ZrSnS and ZrSnSe, respectively. 
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FIG. 11: (Color online) The band structure for (left) 3D and (right) 2D HfSiM with 
M=0, S, Se and Te. The experimental crystal structure of HfSiS is used for HfSiO. For 
others the experimental crystal structure is used. 
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FIG. 12: (Color online) The band structure for (left) 3D and (right) 2D HfGeM with 
M=0, S, Se and Te. The experimental crystal structure of HfSiS is used for HfGeO. For 
others the experimental crystal structure is used. 


20 






















































































HfSnO 



HfSnSe 






fl 




HfSnS 



HfSnTe 



FIG. 13: (Color online) The band structure for (left) 3D and (right) 2D HfSnM with 
M=0, S, Se and Te. The experimental crystal structure of HfSiS is used for HfSnO and 
HfSnS. That of HfGeSe and HfGeTe is used for HfSnSe and HfSnTe, respectively. 
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TABLE I: Binding energy of ZiHM {H=Si, Ge, Sn and M=0, S, Se and Te) with GGA 


Bulk Energy (eV) Single layer Energy (eV) Binding Energy (eV) 

ZrSiO 

-46.130021 

-44.963590 

1.166431 

ZrSiS 

-43.368287 

-42.910675 

0.457612 

ZrSiSe 

-41.300060 

-41.087021 

0.213039 

ZrSiTe 

-37.549321 

-38.490033 

0.059288 

ZrGeO(ZrSiO) 

-43.800281 

-42.645101 

1.155180 

ZrGeS 

-41.324696 

-40.720312 

0.604384 

ZrGeSe 

-39.401731 

-39.076459 

0.325272 

ZrGeTe 

-37.200315 

-37.028521 

0.171794 

ZrSnO(ZrSiO) 

-37.118496 

-36.040903 

1.077593 

ZrSnS(ZrGeS) 

-36.269864 

-35.743714 

0.526150 

ZrSnSe(ZrGeSe) 

-35.139562 

-34.880462 

0.259100 

ZrSnTe 

-35.395074 

-34.708958 

0.686116 


TABLE 11; Binding energy of HfiLM (iL=Si, Ge, Sn and M=0, S, Se and Te) with GGA 



Bulk Energy (eV) Single layer Energy 

(eV) Binding Energy (eV) 

HfSiO(HfSiS) 

-48.066551 

-46.724947 

1.341604 

HfSiS 

-45.837118 

-45.424545 

0.412576 

HfSiSe 

-43.643661 

-43.434275 

0.209386 

HfSiTe 

-41.276397 

-41.249520 

0.026877 

HfGeO(HfSiS) 

-45.697729 

-44.372621 

1.325108 

HfGeS 

-43.692550 

-43.000363 

0.692187 

HfGeSe 

-41.681841 

-41.251401 

0.430440 

HfGeTe 

-39.368634 

-39.127598 

0.241036 

HfSnO(HfSiS) 

-39.050038 

-37.800096 

1.249942 

HfSnS(HfSiS) 

-36.924244 

-36.620387 

0.303857 

HfSnSe(HfGeSe) 

-37.274366 

-36.910855 

0.363511 

HfSnTe(HfGeTe) 

-35.721706 

-35.527372 

0.194334 
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TABLE III; Binding energy of ZrSiM (M=0, S, Se and Te) with LDA 


Bulk Energy (eV) Single layer Energy (eV) Binding Energy (eV) 


ZrSiO 

-50.057196 

-48.554224 

1.502862 

ZrSiS 

-47.962043 

-46.940747 

1.021296 

ZrSiSe 

-45.950533 

-45.203211 

0.747322 

ZrSiTe 

-43.032484 

-42.695313 

0.337171 
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